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ABSTRACT 

We describe a new formalism to fit the parameters a and (3 that are used in the SALT2 model 
to determine the standard magnitudes of Type la supernovae. The new formalism describes the 
intrinsic scatter in Type la supernovae by a covariance matrix in place of the single parameter 
normally used. We have applied this formalism to the Sloan Digital Sky Survey Supernova 
Survey (SDSS-II) data and conclude that the data are best described by a = 0.1351^7 and 
(3 = 3.19^q 24 , where the error is dominated by the uncertainty in the form of the intrinsic 
scatter matrix. Our result depends on the introduction of a more general form for the intrinsic 
scatter of the distance moduli of Type la supernovae than is conventional, resulting in a larger 
value of /3 and a larger uncertainty than the conventional approach. Although this analysis results 
in a larger value of f3 and a larger error, the SDSS data differ (at a 98% confidence level) with 
(3 = 4.1, the value expected for extinction by the type of dust found in the Milky Way. We 
have modeled the distribution of supernovae la in terms of their color and conclude that there is 
strong evidence that variation in color is a significant contributor to the scatter of supernovae la 
around their standard candle magnitude. 
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1. Introduction 

Type la supernovae (SN) have been extensively 
used as "standard candles" to measure astronom- 
ical distances. However, they are often called 
"standardizable candles" to emphasize the fact 
that the light curve shape and color can be used 
to correct the apparent magnitudes, resulting in a 
more accurate measurement of the distance. The 
correlation between light c urve shape and magni- 
tude was demonstrated by I Phillips et al.l (J1993J) , 
who parameterized the light curve shape with the 
parameter Amis , the decay in B-band magnitude 
15 days after the peak. Different parameteriza- 
tions of the light curve shape parameter have been 
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used: A for MLCS2k2 (I Jha et all 120071) . x x for 
SALT2 (ICuv et alfeOOTl feoiOlV and s B for SiFTO 



(jConlev et al.l 120081 ). but the parameter is essen- 



tially determined empirically to produce the tight- 
est correlation between light curve shape and dis- 
tance. This work will use the SupcrNova ANAlysis 
(SNANA) software package (JKessler et al.ll2009ah 
implementation of the SALT2 model that uses x\ 
as the light curve shape parameter. 

Type la SN color depends on the light curve 
shape parameter, but SN also exhibit variations 
in color for fixed light curve shape. The source 
of color variation is not well understood, but it is 
widely suspected that variations in the explosion 
process, host galaxy dust, and the circumstellar 
environment all play a role. The difference in mea- 
sured color and the mean color for a given x\ is 
denoted by the parameter c. Regardless of the pro- 
cess that produces different colors, it may be ex- 
pected that the apparent i?-band magnitude (tub) 
will show some dependence on color, and that an 
appropriate color correction should improve the 
precision of SN as standard candles. These con- 
siderations have led to the widespread practice of 
fitting SN light curves with a one parameter fam- 
ily of templates plus a color law to describe the 
color. 

The determination of color is straight-forward 
given multi-band photometric measurements and 
a color law, and the result is traditionally given as 
E(B — V). However, the amount of correction to 
apply to a given SN has been less clear. The cor- 
rection has been often expressed as an extinction 
correction, which is related to the color through 
the parameter Ry — Ay/E(B — V), where Ay 
is the extinction in ]/-band. Many efforts to 
determine this correction have relied on the pa- 
rameterization ofCardelli, Clayton, and Mathis 
(|Cardelli et al.lll989h . which we will call the CCM 
model, where the parameter Ry changes both the 
shape of the extinction curve as a function of wave- 
length and the amount of extinction for a given 
color excess E(B — V). The simplest approach 
(jjha et al.l 120071 ) is to fix Ry = 3.1, the average 
value determined from the studies of extinction 
in the Milky Way (|Sneden et al.lll978l ). However, 
many studies of SN light curve data have indicated 
that smaller values of Ry provide a better descrip- 
tion of the SN data. In the CCM model the single 
parameter Ry can be measured by estimating Ay 



or by comparing the relative extinction in different 
filters (by considering E(B - V) - E(V - R), for 
example). 

The SALT2 model parameterizes the SN rest- 
frame spectrum S(X, p) as a function of wavelength 
(A) and phase (p) according to 

S(X,p) = xo [S (X,p)+x 1 S 1 (X,p)} e [- cCL Wl, (1) 

where the fitted parameters xo, %i, and c are the 
overall scale, the light curve shape parameter, and 
the color respectively. The functions So and Si are 
fixed by a training sample of SN light curves and 
spectra. The SALT2 light curve fit does not rely 
on the CCM parameterization, but uses a color law 
CL(X) determined empirically from SN data. The 
SALT2 parameters are determined by fitting the 
measured light curve data for each SN to the syn- 
thetic magnitudes calculated from the redshifted 
spectrum given by Equation ([TJ . 

The light curve parameters xo, x%, and c spec- 
ify a model SN spectrum, but an additional step 
is needed to use the light curve parameters to cal- 
culate a standard SN magnitude. One can calcu- 
late the apparent magnitude tub , where tub is the 
apparent _B-band synthesized from the rest-frame 
spectrum given by Equation (Q]), and a standard- 
ized SN brightness m s B is computed with parame- 
ters a and /3: 



= rriB + olx\ — f3c. 



(2) 



iGuv et alJ (|2007l ) describe a technique that si- 
multaneously determines the cosmological param- 
eters and the coefficients (a and /3) that mini- 
mize the residuals on the Hubble diagram. The 
SALT2 approach to correction for color is formally 
equivalent to the description in terms of Ry and 
E(B — V), except for the interpretation that the 
value of (3 may result from some effect other than 
extinction by host galaxy dust. Regardless of the 
interpretation, /3, which gives the magnitude cor- 
rection for i3-band, and Ry, which gives the mag- 
nitude correction for V-band, may be compared 
via the relation /3 = Rb = Ry + 1. 

In this work, we present a new method for fit- 
ting the parameters a and (3 from the data. The 
method has been implemented in the program 
SALT2mu, which has been added to the SNANA 
software package. The motivations for a new pro- 
gram are: 



• The fit for a and (3 is decoupled from the cos- 
mology. If the mean values of x\ and c are 
a function of redshift (a magnitude limited 
survey will have such selection effects), the 
average residual depends on both the cos- 
mology and the a and /3 parameters. We 
prefer to determine the a and j3 solely by 
minimizing the rms of the distribution and 
to let the mean of the distribution determine 
the cosmology. 

• A new description of intrinsic scatter is used 
as discussed below. 

• The SALT2mu output can be used inter- 
changeably with different cosmology fitters 
since SALT2mu produces distance moduli 
separately from the cosmology fit. This fea- 
ture allows a direct comparison of SN dis- 
tance moduli without reference to a partic- 
ular cosmology. 

• The sample used to determine a and /3 can 
be different from the sample used to deter- 
mine the cosmology. 

A significant additional benefit to writing a new 
program is that SALT2mu reads files from the 
SALT2 light curve fits and outputs files that can 
be read by any of the SNANA cosmology fitters. 
The fitting model resides in a single function mak- 
ing it easy to modify the \ 2 function by adding 
new parameters, for example. 

The inclusion of a more general description of 
the intrinsic scatter of the fitted SALT2 parame- 
ters is a major feature, as is the determination of a 
and /3 in a way that is independent of the cosmol- 
ogy. The details of the incorporation of the intrin- 
sic scatter are given in the mathematical descrip- 
tion ($2]) below, and most of the paper is devoted 
to the issues associated with the more general de- 
scription of the intrinsic scatter. The program is 
then applied to the Sloan Digital Sky Survey su- 
pernova (SDSS-II SN) data, which is used to il- 
lustrate some of the features of the program and 
issues with its use. 



The SDSS-II SN data ( Frieman et al.1 2008) 



were obtained during 2004-20 07 as part of ex - 
tension of the origi nal SPSS (lYork et al.l l2000l). 
The SDSS t elescope (iGunn et al.ll2006l ) and imag- 
ing camera (Gunn et al. 1998) produce photomet- 
ric measurements in each of the 5 SDSS filters 



(JFukugita et al.lll996l) spanning the range of 350 
to 1000 nm. The most useful filters for observing 
SDSS SN, however, are g, r, and i because the 
SN are too faint in u and z to be well measured 
by the SDSS SN survey except at low redshifts. 
Further details of the survey, data processing and 
selections are given in $3] 

2. Mathematical Description 

The standardized 10 pc supernova magnitude 
is determined by the expression 



M, 



(B) _ 



tub — n(z) + a(z)x\ — /3(z)c, (3) 



where /i(z) is the distance modulus, and Mq is 
the standardized SN magnitude for B-b&nd. We 
have generalized the unknown parameters a and /3 
to functions of redshift to accommodate possible 
effects of evolution of supernovae or their environ- 
ment. These functions are not known from first 
principles, but must be determined from the data. 
We use a nearly equivalent formulation 



M (z) = m x - /i(z) + a(z)xi - P(z)c, (4) 



where 



-2.51og 10 (a;o) 



(5) 

The definition shown in Equation ([5]) has the ad- 
vantage that the magnitude m x does not depend 
on the shape of the i?-band filter. It differs from 
the _B-band standard magnitude by about 10 mag- 
nitudes. The functions Mq(z), a(z) and /3(z) are 
to be found by the minimization of the expression 
for x 2 

N 

X 2 = ^ [m x „ - n(z„) + a{z n )x ln - f3(z n )c n 

-M (z n )] 2 /(a 2 n +a int ) 2 (6) 

where m xn , x\ n and c„ are the fitted SALT2 pa- 
rameters for the n th SN, and N is the number of 
SN. The quantity a n is the error on the quantity 
in the numerator due to measurement uncertainty 
for the n th SN. The quantity eri nt is discussed 
below and its form is displayed in Equation (0. 
The function M (z) is, by construction, degener- 
ate with cosmological information contained in the 
distance modulus fi{z). The introduction of a nui- 
sance function Mq(z) allows our fit to minimize the 



scatter in the Hubble diagram while being insensi- 
tive to the cosmology. In order to reduce this func- 
tion to a few parameters we assume that Mq(z) is 
an arbitrary constant in each of several redshift 
bins centered at z — Zb, and we use the exact ex- 
pression for the distance modulus fj,(z) for some 
assumed cosmology. Effectively we are using an 
assumed cosmology to extrapolate the residual to 
the redshift at the center of the bin and minimiz- 
ing the Hubble diagram scatter at those reference 
points. Since our method introduces several new 
parameters, the statistical accuracy of the fit is 
reduced by the additional parameters. And since 
our method relies on the minimization of the scat- 
ter of distance moduli in a single bin, it requires 
a bin population of at least two SN — preferably 
many more. These considerations are of no conse- 
quence for the SDSS sample, but may be of some 
concern in applying our method to samples with 
sparse data or unconventional cosmologies. We 
show in <g] that the fit results are insensitive to 
the choice of cosmological parameters if the num- 
ber of redshift bins is sufficiently large. 

Although this formalism allows fitting an arbi- 
trary redshift dependence! for a and /3, we will 
assume that a and j3 are constants, independent 
of redshift , in this paper. Recently, a number 
of papers (Kelly et al.ll2010t iLampeitl et alJl2010t 
Sullivan et alJl2010l ) have found that the value of 



Mq shows a dependence on host galaxy proper- 
ties. Other papers have found corr elations be- 



tween line velocities in S N spectra ([Wang et al 



2009; Foley & Kascn 2011) and their standardized 
magnitu des and other corre lations with spectral 
features ( Nordin et alJl201lh . Those aspects will 
not be explored in this paper, but it is straightfor- 
ward to generalize Equation Q by adding a new 
term "/y n to the numerator of Equation ([6]) , where 
y n is some measured property of the n SN and 7 
is a parameter to be determined by the minimiza- 
tion. In addition, the measurement uncertainty 
and intrinsic scatter in y n should be incorporated 
into <r„ and crj nt , respectively, following the for- 
malism described below. 

When the SALT2 light curve parameters are 
used according to Equation (j4j to construct a SN 
Hubble diagram, the dispersion exceeds that ex- 



1 SALT2mu includes the possibility of fitting a linear depen- 
dence for a and 0. 



pected from the light curve fit parameter errors. 
A simple way to parameterize the Hubble diagram 
scatter is to assert that the distribution of light 
curve parameters can be described by an arbitrary 
population in the plane described by Equation (UJ) 
plus a Gaussian spread that is described by a co- 
variance matrix, which we will refer to as the in- 
trinsic scatter. The quanitity CTi n t in the denomi- 
nator of Equation ((6} describes the expected error 
in the numerator due to the intrinsic scatter. We 
further assume that the covariance matrix is con- 
stant, independent of the SN light curve parame- 
ters. These considerations result in an expression 
for the intrinsic scatter 

af nt = S 00 + a(z) 2 £ n + /3(z) 2 £ cc + 2a(z)£ i 
- 2/3(z)E 0c - 2a{z)/3(z)Z lc , (7) 

where £ is the intrinsic covariance matrix of the 
parameters and £00, ^11 1 an d ^cc are the diagonal 
elements that correspond to the SALT2 parame- 
ters m x , x\, and c. Our introduction of a covari- 
ance matrix to describe the intrinsic scatter differs 
from the conventional approach where the intrinsic 
scatter is assumed to be a single number. Equa- 
tion ([7]) shows that the conventional approach is 
equivalent to our model if £00 ^ s the only non-zero 
term, leading to the conclusion that the conven- 
tional approach applies the intrinsic scatter to m x 
but not to X\ or c. We introduce the more gen- 
eral approach because the additional terms have 
not been excluded on either experimental or the- 
oretical grounds and, in fact, there is good rea- 
son to believe that there is a significant scatter 
that can be attributed to the SALT2 color param- 
eter. The measurement error term (a 2 ) in Equa- 
tion (|6]) is dominated by photometric errors and 
is calculated by the SALT2 light curve fit for each 
SN and is given by a covariance matrix that is of 
the same form as the intrinsic scatter covariance 
displayed in Equation ([7]). The denominator in 
Equation (|6]) is therefore the sum of a constant in- 
trinsic scatter covariance matrix, which we have 
introduced above, and the measurement covari- 
ance, which has traditionally been included. The 
fact that the measurement error and the intrin- 
sic error are both of the form shown in Equation 
([7]) means that each of the SALT2 parameters has 
both an intrinsic scatter and a measurement error 
that are taken into account in computing y 2 . The 
formulation of the minimization problem defined 



by Equations (|6]) and Q is an ex tension of the 
formula given in IPress et al.l ( 20071 ) for the case of 
fitting a line to two variables with known errors. 
A more extensive treatment of the linear regres- 
sion problem from an astron omical persp ective has 
been given, for example, by iKellvl ( 2007 ). 

After the parameters a and (3 are determined 
by the minimization of x 2 the distance moduli arc 
calculated according to 

Mn = m xn + a{z n )x ln - /3(z n )c n - M (z b ) (8) 

These distance moduli can then be used for input 
to cosmological fits. 

3. Application to the SDSS data 

We now apply this method to the SDSS 4-year 
sample of 529 spectroscopically confirmed SN. The 
SDSS SN survey is a so-called "rolling search," 
where a portion of the sky is repeatedly scanned 
to discover new SN and to measure the light curves 
of ones already discovered. The survey measured 
an equatorial stripe about 2.5 degrees wide in dec- 
lination between a right ascension of 20 h and 04 
h. Full coverage of the stripe can be obtained in 
two nights, but the average cadence was about 
four nights because of inclement weather and in- 
terference from moonlight. The data were ob- 
tained in the Fall of 2005, 2006, and 2007 and a 
smaller amount of data was taken in a test run in 
2004. The survey is sensitive to SN beyond a red- 
shift of 0.4, but beyond a redshift of 0.2 the com- 
pleteness and the ability to obtain accurate pho- 
tometry deteriorates. The SDSS camera images 
were processed by the SDSS imaging software and 
SN were identified via a frame subtraction tech- 
ni que. More deta i ls and references can be found 
in iFrieman et al.l ( 20081 ). The candidate selec- 



using techniques s imilar to those use d for the 
first year sample (jZheng et al.l 12008) 



tion and s pectroscopi c ident ification has been de- 
scribed bv lSako et al.l (|2008l) . The supernova pho- 



tometry used in this analysis wa s obtained using 



the sc ene modeling technique of iHoltzman et al 



(|2008l ). where the sequence of SN observations is 
modeled as a variable point source and a host 
galaxy background that is constant in time. 

We select a sample of well-measured light 
cu rves for this analysis b ased on the criteria used 
by iKessler et al.l (|2009bl ) . We require 



• Redshift in the range 0.02 < z < 0.42 

• A minimum number of measurements with- 
out any requirement on the measurement er- 
ror and a minimum number of well-measured 
points, where well- measured means an error 
of 20% or less {i.e., S/N > 5) 

— At least 5 well-measured points 

— At least 2 well-measured points in dif- 
ferent filters 

— At least 1 point measured before peak 
light in B-band 

— At least 1 point measured later than 9.5 
days in the rest frame after peak light 
in B-band 

• An acceptable fit to the SALT2 model (con- 
fidence level greater than 0.001) 

These selections result in a total sample of 343 
SNe with negligible contamination from sources 
that are not Type la SN. The redshift range spans 
0.036 to 0.419, with a mean redshift of 0.22. The 
SDSS sample thus defined does not have uniform 
selection criteria as a function of redshift and, 
therefore, there is a significant change in the pop- 
ulation as a function of redshift. In particular, 
the spectroscopically-selected sample of SDSS SN 
is significantly biased towards bluer, brighter SN 
as the redshift increases; the sele ction b ias is dis- 

(l2009bl) . To the 



cussed in detail in Kessler et al. 



Spectroscopically confirmed as a Type la SN 



extent that the data are described by the model of 
Equation Q with an accurate parameterization of 
a{z) and (3(z) and the SALT2 fit parameters are 
unbiased, we expect to be insensitive to having a 
biased population. 

An application of SALT2mu to the SDSS data 
using four redshift bins yields the results shown 
in Table [U which shows the fit x 2 for 337 degrees 
of freedom (343 SN), the fitted values of a and 
/3, and their errors a a and a p. In section §4 we 
show that four bins is adequate to decouple the 
SDSS data from the determination of the cosmo- 
logical parameters. The first row in Table [T] shows 
the fit when the intrinsic scatter is assumed to be 
zero resulting in \ 2 — 1305, reproducing the well- 
known result that the scatter in the Hubble dia- 



gram is larger than what is expected from the mea- 
surement error alone. Subsequent rows illustrate 
how the fitted parameters change with different 
assumptions on the nature of the intrinsic scatter. 
Each of the fit variables is assumed, in turn, to be 
totally responsible for the intrinsic scatter with the 
magnitude of the scatter being adjusted to achieve 
X 2 /d°f ~ 1. In Tabled] we have introduced new 
variables to indicate the size of the intrinsic scat- 
ter: Co = n/Sooj o"i = VEu, and cr c = \/E cc . 
From Table Q] we see that the different assump- 
tions on the nature of the intrinsic scatter result in 
rather different estimates of the a and /3 parame- 
ters compared to the reported statistical error. An 
assumed scatter in m x (<r ^ 0) decreases both a 
and /? compared to the case where the intrinsic 
scatter is taken to be zero. An assumed scatter in 
x\ (ci 7^ 0) increases a by 15.7er while decreas- 
ing j3 by 1.4<7 relative to the fit where the intrinsic 
scatter is assumed to be entirely due to m x . The 
opposite effect is seen when the scatter is assumed 
to be due to c (<r c ^ 0): /3 increases by 7.1cr and a 
decreases by 3.1er. Even more variation could be 
exhibited by a more general intrinsic scatter co- 
variance matrix. However, the variations in a and 
/3 are not without limit. The values are all the 
same sign and approximately the same magnitude. 
Roughly speaking, we are changing our estimate 
of how much to correct the observed slope for the 
flattening caused by the scatter described by the 
intrinsic covariance matrix (Equation ([7]l). 

If we are only interested in SN as distance in- 
dicators, a and /3 are nuisance parameters of no 
particular interest. The average SN brightness at a 
given redshift will be given by the mean of Equa- 
tion @: if the mean values of x± and c do not 
evolve with redshift, then different values of a and 
j3 will result in a redshift-independent offset that is 
not significant in determining the cosmology. The 
mean values of x\ and c for the SDSS data are 
plotted in Figure Q] showing, however, a substan- 
tial increase in x\ and a decrease in c with red- 
shift. These trends are primarily caused by the 
spectroscopic-selection bias towards brighter SN 
la. The different values listed in Table Q] would 
result in differences of ~ 0.1 in the average val- 
ues of fx(z), a rather large error for cosmological 
studies. In practice, the effect can be mitigated by 
limiting the sample to lower redshifts, where the 
selection efficiency is higher, correcting for biases, 



and combining high and low redshift experiments 
that have different efficiencies. Nonetheless, the 
possibility of an uncorrected bias is a concern. 

Our results can be c ompared to the results of 
Amanullah et al.l ( 20101) who combined data from 
several sources to obtain a value of/3 = 2.51±0.07 
with an an alysis most simila r to our case with 
er = 0.14. iGuv et al.1 (l2007h found /3 = 1.77 ± 
0.16 for the fi rst yea r data from the SNLS and 
Kessler et all (l2009bl) found /3 = 2.65 ± 0.08 for 



Amanullah et al 



the fi rst year of SDSS data 
(|2010l ) also quote numbers for their four largest 
samples: /3 = 2.38 ± 0.14 (SDSS), 2.73 ± 0.13 
(CfA), 2.50 ± 0.17 (E SSENCE), and 1.72 ± 0.18 
(SNLS). More recently IConlev et al.l (J201ll ) found 
(3 = 3.18 ± 0.10 for the SNLS 3-year sample us- 
ing a chi-squared minimization technique. The 
study of lNobili and Go obar (2008) found the vari- 
ations in SN colors were most consistent with a 
CCM param eterization when Ry = 1 .75 ± 0.27. 
However, the lNobili and Goobarl ( 2008 ) study dif- 
fers significantly from this analysis in that their 
value of Ry is based on their measurement of the 
wavelength dependence of t he reddening using the 



Kessler et al. 



(|2009bl) 
2.18± 



Cardelli parameterization, 
also investigated SN color and found Ry 
QAA(stat) ± 0A8(syst), but that result relies on 
modeling the color changes due to the magnitude 
limited sample. All these results are broadly con- 
sistent in finding f3 < 4.1, but they exhibit more 
scatter than one would expect from the quoted er- 
rors. 

4. Tests with simulated SN samples 

We will return to the SDSS SN data in g5j but 
it is useful to first examine the performance of 
SALT2mu on simulated SN data. We first deter- 
mine the number of data bins that are required to 
decouple our fits from the cosmology over a wide 
range of parameters and conclude that 4 bins are 
adequate for the SDSS data. We then verify the 
formalism with a simplified simulation and show 
that we recover the input values of a and j3 with- 
out bias. We next look at a realistic simulation 
and show that the fitted a and /3 dependence on 
the intrinsic covariance matrix £ follow the trends 
seen in the data. We compute the intrinsic scat- 
ter covariance matrix from the simulation (which 
requires knowing the simulation parameters), and 



Table 1 
SALT2mu fits to the SDSS Data 





Intrinsic Error 




x 2 


a 


c a 


P 


C/3 


CO 


= Ci = 


Cc = 


= 




1305 


0.1891 


0.0064 


2.820 


0.050 


Co 


= 0.14 


Cl = 


c„ 


= 


350 


0.1465 


0.0098 


2.650 


0.078 


CI 


= 0.70 


CO = 


Co 


= 


330 


0.3002 


0.0200 


2.405 


0.116 


Gc 


= 0.05 


CO = 


Cl 


= 


333 


0.1330 


0.0110 


3.208 


0.102 



Note. — The SALT2mu program is applied to the 343 selected SDSS 
SN. The values of a and /? and their calculated errors {a a and (Jp) are 
shown. There are 6 fit parameters (a, (3, and the 4 Mq values for the 4 
redshift bins) so the number of degrees of freedom is 337. The variable 
Co, Ci, and a c give the assumed rms scatter in the SALT2 parameters 
m x , x±, and c, respectively. The first row in the table shows the fit 
assuming an intrinsic scatter of zero; the following 3 rows assume that 
the intrinsic scatter results, in turn, from a single SALT2 fit parameter. 




Redshift 



Fig. 1. — The average value of (a) SALT2 shape 
(xi) and (b) color (c) as a function of redshift for 
the SDSS sample of 343 SN. The general trends of 
the data can be explained by selection effects. 



show that the input values of a and /3 are recov- 
ered when the correct intrinsic scatter covariance 
matrix is used. 

We generate the simulated SNe by using 
SNANA to generate SN light curve data according 
to the SALT2 model. The SN are generated with 
a redshift distribution appropriate for the SDSS 
survey and with a — 0.18 and /3 — 3.4. Realistic 
cadence, measurement errors, and selection effects 
are applied to the simulated SN, and they are pro- 
cessed with the same light curve fitting program 
and selection criteria as the SDSS SN, resulting in 
SALT2 fit parameters for the simulated SN. In ad- 
dition to the photometric measurement errors, the 
simulation applies color smearing to model the in- 
trinsic scatter. Color smearing consists of random 
magnitude offsets that are generated for each SN 
and observed filter and added to each light curve 
point measured in that filter. The magnitude of 
the smearing is adjusted to 0.08 in each observed 
filter to give approximately the same Hubble dia- 
gram scatter as is seen in the real data. 

The first test is a relatively trivial one to verify 
that the fit results for a and j3 are, indeed, in- 
dependent of cosmology. We expect that a large 
number of parameters Mo(zb), will decouple the 
determination of a and (3 from the cosmological 
parameters, but wish to minimize the number of 
nuisance parameters that are introduced. We il- 
lustrate the fit dependence on the number of red- 
shift bins in Figure [5] using the simulation of the 
SDSS dataset, which has over ten times the num- 
ber of SN (4417) as the SDSS data. Figure |2] 



shows the fit results for a and /3 assuming three 
different cosmologies: VLm = 0.27, 51a = 0.73, 
w = -1 (red), M = 0.27, J1a = 0.0 (green), and 
n M = 0.27, 0\ = 0.73, w = (blue). The differ- 
ence in the fitted parameters is noticeable when 
only 1 bin is used, but becomes negligible as the 
number of bins becomes greater than about 5. We 
conclude that using 4 redshift bins is adequate for 
the SDSS data, given the lower statistical precision 
of the data and the established constraints on the 
cosmological parameters. All results in this paper 
are based on dividing the SDSS data (real or simu- 
lated) into 4 equal sized bins spanning the redshift 
range 0.02 < z < 0.42. 

While the fit results in Figure [5] show that we 
have successfully decoupled the fit for a and /? 
from the cosmology, the dependence on the as- 
sumptions about the intrinsic scatter matrix re- 
mains. For the fit results shown in Figure [5] we 
took (To = 0.14 and <j\ = o c = 0, and this choice 
results in fitted values of a = 0.184 and j3 — 2.86 
which differ from those used to generate the sim- 
ulated SNe, namely a = 0.18 and j3 = 3.4. The 
discrepancy is not surprising since our simulation 
includes significant color smearing that we have 
not taken into account. We will determine more 
accurate parameters for the intrinsic scatter co- 
variance below, but first we check the formalism 
with a simplified simulation. We use the same 
simulated SN but replace the simulated values of 
m x , x\, and c with the exact values that were used 
to generate the light curves but smeared accord- 
ing to the 3x3 error matrix including any as- 
sumed intrinsic scatter. This procedure tests the 
self-consistency of the SALT2mu model and elimi- 
nates concerns about biases in the SALT2 param- 
eter fitting and inaccuracies in the calculation of 
the parameter errors. The results of the fits are 
shown in Table [2] with the same four assumptions 
about the intrinsic scatter that were used in Table 
[TJ We have also tested the formalism when all the 
elements of the intrinsic scatter covariance matrix 
are non-zero; the last line in Table [2] is an exam- 
ple. For each case we recover the input values that 
are consistent within the statistical error of about 
1%. The x 2 /dof consistent with 1. We do not ob- 
serve any variation in the parameter values since 
the parameters of each SN are generated according 
to a Gaussian error distribution that is the com- 
bination of the measurement error and whatever 
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Fig. 2. — The dependence of the fitted parame- 
ters a (a) and (3 (b) is displayed as a function of 
the number of redshift bins for three assumptions 
about cosmological parameters and a particular 
choice of the intrinsic scatter matrix (crn = 0.14, 
o\ = cr c = 0). As the number of bins increases, 
the fits converge to a common value a = 0.184 and 
f3 = 2.86 shown as dotted lines on the respective 
figures. 



Table 2 
salt2mu fits to the simplified, simulated sdss data 



Intrinsic Scatter 


x 2 


a 


CT a 


P 


a H 


(To — <J\ — <7 C — 


4445 


0.1802 


0.0009 


3.391 


0.010 


CTo = 0.14 


4487 


0.1774 


0.0019 


3.395 


0.024 


eri = 0.70 


4319 


0.1815 


0.0018 


3.423 


0.023 


a c = 0.05 


4277 


0.1818 


0.0022 


3.379 


0.028 


Full covariancc a 


4564 


0.1811 


0.0015 


3.356 


0.019 



a Thc full covariancc matrix used for this example was ctq — 0.05, 
rxi = 0.35, er c = 0.025, £oi = -0.2, { c = 0.5, and £ lc = 0.4, where 
the £ are the correlation coefficients defined in the text following 
Equation ( f9l . 



Note. — In the SALT2mu simplified simulation, the SALT2 fit 
parameters (m^,, x\, and c) arc replaced with the true parame- 
ters plus an random error determined by the combination of the 
measurement and intrinsic scatter. The number of simulated SN 
is 4417 and the number of degrees of freedom for this fit is 4411. 



intrinsic scatter was assumed. 

Using the parameters of the SALT2 fits to the 
simulated data provides a more complicated test 
since the color smearing model is different than ap- 
plying an intrinsic scatter to the simulated SALT2 
parameters used to generate the light curves. We 
apply SALT2mu to the simulated SN using the 
same procedure used for the SDSS data, and dis- 
play the results in Table El The first line has 
the intrinsic scatter set to zero, and the succeed- 
ing lines have the intrinsic scatter that results in 
X 2 /dof ~ 1. The general trends shown in Table [3] 
are very similar to those seen in Table [TJ assum- 
ing an intrinsic scatter in x\ and c increases the 
fitted value of a and /3, respectively. The magni- 
tudes of the intrinsic scatter required to achieve 
X 2 /dof sw 1 are similar to those required by the 
SDSS data. But we can also compare the re- 
sults with the "true" values for the simulated data. 
While the assumption that the intrinsic scatter is 
entirely due to color errors (last line in Table [3]) 
produces values of a and f3 that are close to those 
used to generate the simulation (a differs by 2.6a), 
it is not a very accurate representation of the in- 
trinsic scattering covariance matrix that describes 
the simulated SN as shown below. 

In principle, we know the smearing formulae 
applied in the simulation and can calculate the 
intrinsic scatter covariance matrix E. But the cal- 
culation of the intrinsic smearing matrix is compli- 
cated because the error model used in the simula- 



tion smears individual light curve points in a way 
that is not easily expressed in terms of an intrin- 
sic smearing matrix in m x , Xi, and c. However, it 
is relatively straightforward to calculate the error 
a posteriori. The full error matrix is determined 
by the variance of the simulated measurements of 
m x , x\, and c relative to the true values (which, 
of course, are known for the simulated data). The 
light curve fitter calculates covariance matrix of 
measurement errors, which, when added to the in- 
trinsic scatter covariance matrix, should yield the 
total error matrix. The intrinsic scatter covari- 
ance matrix can therefore be calculated for the 
simulated data by subtracting the observed vari- 
ance from the prediction for the measurement er- 
ror. The total error covariance matrix calculated 
from the variance of the reconstructed SALT2 pa- 
rameters from the simulated ones is 



ao 


= 0.0781, 


G>i 


= 0.253, 


£oc : 


o-i 


= 0.351, 


£ic 


= 0.183, 


and 


O r 


= 0.0634. 









0.612, 



(9) 



where the various £ are the correlation coefficients 
(£oi = Soi/((Tocri), £oc = E 0c /(ctoO'c), and £i c = 
Ei c /(ericr c )). Using the same notation to represent 
the elements of the measurement error matrix we 
find 



a = 0.0365, 


£oi 


= 0.581, 


£oc : 


o-i = 0.347, 


£ic 


= 0.307, 


and 


a c = 0.0286. 









0.694, 



(10) 



Table 3 
salt2mu fits to the simulated sdss data 



Intrinsic Error 



X 



a 



CT/3 



<7q — ai — a c — 
cto = 0.14 
cti = 0.65 
ac = 0.045 



19500 
4340 
4375 
4343 



0.2143 
0.1846 
0.2482 
0.1855 



0.0011 
0.0019 
0.0025 
0.0021 



3.104 0.011 

2.868 0.021 

2.882 0.026 

3.370 0.025 



Note. — The first row in the table shows the fit assuming an in- 
trinsic scatter of zero; the following 3 rows assume that the intrinsic 
scatter results, in turn, from a single SALT2 fit parameter. The 
number of simulated SN is 4417 and the number of degrees of free- 
dom for this fit is 4411. 



co 


= 0.0690, 


&i 


= -0.119, 


£0c = 


Cl 


= 0.0506 


6c 


= 0.362, 


and 


<7 r 


= 0.0566. 









Subtracting (fl~0|) from ((9j gives an excess error, 
which we attribute to the intrinsic scatter. The 
result is 



0.590, 



(11) 
We are now in position to apply the intrinsic 
scatter from Equation (fTTj) above to the fit from 
SALT2mu. The result is shown in the first row in 
Table QJ The result is consistent with the gener- 
ated values (a = 0.18 and /3 = 3.4). The small 
discrepancy in j3 (2.6er) not a serious concern be- 
cause it is quite small compared to the accuracy of 
the SDSS dataset. However, we suspect that un- 
certainties in the calculation of the measurement 
errors, including treatment of non-linearities and 
a consistent handling of model errors, may con- 
tribute a systematic bias at this level of precision. 
We have performed additional fits with other in- 
trinsic scatter matrices. The second row is a fit 
with all the non-diagonal elements of the intrin- 
sic scatter matrix set to zero. The fitted param- 
eters are only slightly changed although there is 
a decrease in \ 2 ■ The next row shows a fit that 
has been selected to give a larger change in \ 2 by 
tuning the correlations for maximum effect. The 
last row in Table |4] displays a fit, which may be 
counter-intuitive: the error in x\ is decreased, but 
the x 2 increases. These results illustrate how com- 
plicated the fit behavior can be for the five free pa- 
rameters that live in the intrinsic scatter matrix. 



5. Determination of the intrinsic scatter 
matrix 

We do not have a rigorous approach to estimat- 
ing the intrinsic scatter matrix for actual SN data. 
Since we don't know the underlying, unsmeared 
distribution or the actual form of the smearing, 
the observed distribution is difficult to interpret. 
The scatter in the Hubble diagram sets the overall 
scale, but more information is needed to determine 
the individual elements of the intrinsic scatter ma- 
trix. 

We can derive some information from the color 
distribution with the help of some pl ausible as- 
sumpt ions, which were first used by Ijha et al. 
()2007n . We assume that the color distribution 



arises from random color smearing that we pa- 
rameterize as a Gaussian of unknown width (a c ) 
and host galaxy extinction that we characterize 
as an exponential distribution of unknown slope 
(r). We do not expect to be sensitive to the exact 
shape assumed for galaxy extinction: the impor- 
tant point is that it is one-sided (extinction can 
only make objects redder, not bluer) and that the 
most probable value of extinction is zero. In addi- 
tion, we do not know the position of the edge of the 
extinction distribution, so we must introduce an 
additional parameter (co) to be determined from 
the data. Assuming that the intrinsic scatter and 
host galaxy extinction are independent, the convo- 
lution of the two processes results in the following 
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Table 4 

salt2mu fits to the simulated sdss data 
with non-zero, off-diagonal elements of the intrinsic scatter matrix 



Intrinsic Er 



P 



As computed 

As computed but £01 — £oc — £ic — 0.0 

As computed but | i = 0.847, £ 0c = -0.913, £i c 

As computed but a\ — 





4161 


0.1832 


0.0022 


3.469 


0.027 




2752 


0.1822 


0.0026 


3.320 


0.031 


0.831 


1758 


0.1820 


0.0031 


3.227 


0.036 




3130 


0.1817 


0.0024 


3.400 


0.030 



Note. — The first row shows the intrinsic scatter matrix computed as described in the text and succeeding 
rows illustrate various changes to the intrinsic scatter matrix. The simulation is the same as was used to 
produce Table [3] 



distribution for c: 



p(c) = A / exp 



c - c n - c 



x exp 



2a 



, , dc' (12) 



In order to minimize the contribution of measure- 
ment errors to the width of the Gaussian com- 
ponent, we include only SN where the estimated 
measurement error for c is less than 0.04, reduc- 
ing our sample of 343 SN to 247. Fitting Equa- 
tion dl2]) to the data, we find a c = 0.061 ± 0.007, 
r = 0.076 ± 0.010, and c = -0.073 ± 0.009. The 
data and resulting fit are shown in Figure |3l As 
can be seen, the functional form of Equation (TT2"j) 
provides a good description of the data. The fit- 
ted value for r is substantially sm aller than the 
value r = 0.334 ± 0.088 quoted by iKessler et al 



(J2009bJ), but our result is an effective r for the 
spectroscopically-selected sample with the addi- 
tional requirement of a well-measured color and 
is therefore biased against large extinctions. Mea- 
surement errors are calculated by the SALT2 light- 
curve fit and contribute, on average, y/Jcrf) = 
0.026 to the width of the Gaussian component, 
resulting in an estimate of 0.055 ± 0.007 for the 
color smearing. We do not expect our values of Cq 
and cr c to be significantly affected by the selection 
bias, and we do not make a correction for the bias. 
This result is comparable to the amount of color 
smearing (a c = 0.05) that was required in Table 
[1] to explain, by itself, the scatter in the Hubble 
diagram. 

Unfortunately, we do not have any useful a pri- 
ori model for the x\ distribution. We can, how- 



n 



fit prediction 
fitted curve 



\ 




SALT2 color (c) 

Fig. 3.— The distribution in SALT2 color of 247 
SDSS SN selected to have a calculated measure- 
ment error er c < 0.04 is shown as a histogram, and 
the dashed curve is the best fit to the data. Indi- 
vidual points on the curve are shown at the center 
of the histogram bins with errors corresponding to 
the expected fluctuations from Poisson statistics. 
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ever, make the somewhat trivial observation that 
the intrinsic scatter cannot be larger than the 
width of the distribution. It is known that pas- 
sive and star forming galaxies have different dis- 
tributions in x\ and, if we assume that the in- 
trinsic scatter is the same for both, we can obtain 
an upper limit for the intrinsic scatter from the 
distribution that excludes passive galaxies. The 
distrib utions in x\ were shown by lLampeitl et al 



(|2010h for the SDSS data, and we have used that 
analysis for the 296 SN this paper has in common 
with that analysis. The resulting distribution is 
shown in Figure |4] along with a Gaussian of width 
1.0. The Gaussian is not a good fit to the tails 
of the distribution, but it is a good match to the 
core and we conclude that an upper limit for the 
intrinsic scatter in x\ is 1.0, a limit that is greater 
than the amount required to explain the scatter in 
the Hubble diagram. 

The simulated data exhibits a significant corre- 
lation between m x and c, and it is interesting to 
ask whether that correlation should be expected 
in the data. In the simulation, the correlation 
arises because the SALT2 model defines i?-band 
to be the reference point for the color and because 
the measured data is largely redward of B-band. 
The assumed independent variation of each mea- 
sured band results in the correlation between peak 
magnitude and color. While the model of color 
smearing is arbitrarily constructed to reproduce 
the scatter in the Hubble diagram, there is some 
evidence to support independent fluctuations in 



each of the filters. As pointed out in (Kcsslcr et al 



2010)), if the redshift is fit from the light curve 
data, the dispersion of the results around the spec- 
troscopically determined redshifts is larger than 
can be accounted for by photometric errors alone, 
but the color smearing model of SNANA gives ap- 
proximately the right dispersion. Based on these 
considerations, we conclude that the correlation 
between m x and c is likely to be positive (as seen 
in the simulations) for the SDSS data. 

While we have been able to make some ob- 
servations on the elements of the intrinsic scat- 
ter matrix, our picture is incomplete and rests 
on assumptions of uncertain validity. An alter- 
native approach is to marginalize over all possi- 
ble intrinsic scatter matices. We implement this 
marginalization by running SALT2mu repeatedly 
with randomly generated intrinsic scatter matri- 



ces. We start with a vector randomly generated 
on the unit sphere subject to the restriction that 
all components of the vector are positive. We then 
multiply the three components of the vector by 
0.2, 1.0, and 0.07 respectively and set <7o, ci, and 
<r c to the scaled vector components. The constants 
are chosen to be about 50% larger than the val- 
ues (see Tabled]) that were required to explain the 
intrinsic scatter. The normalized correlations are 
chosen randomly between —1 and +1. 

A total of 10000 intrinsic scatter matrices were 
generated, but the procedure does not guarantee 
that the matrices are positive. A total of 6380 of 
the matrices were found by SALT2mu to be posi- 
tive definite and were used to determine values for 
a and /3. Each intrinsic scatter matrix is scaled 
by an arbitrary factor to achieve a x 2 /dof = 1. 
The red (dotted) curve in Figure [5] shows the re- 
sults of the fits for (a) a and (b) /?. The median 
value is a = 0.150±g;gl? and f3 = 2.91±q j°, where 
the error is the range that includes 68% of the 
results. Our intrinsic scatter covariance matrix 
defines an intrinsic error in 3 dimensions, but the 
chi-squared is only sensitive to the error in one 
direction. The intrinsic errors in the other two di- 
rections just broaden the distribution of points in 
the plane defined by Equation ((8} ■ If the intrinsic 
errors are large enough to explain, by themselves, 
the observed population, the slopes in Equation 
([8]) become undefined and produce the outliers 
that are observed in Figure [SJ The data show a 
secondary peak near j3 = 2.5, but this should not 
be interpreted as evidence for 2 populations with 
different slopes. The shape of the plot in Figure 
5, of course, depends on the distribution of covari- 
ance matrices, which we chose to be uniform. The 
lower values of j3 result primarily from intrinsic 
scatter matrices with a small a c . The fit becomes 
insensitive to a c as a c — > and the peak at ft = 2.5 
is an artifact of our choice of distribution for the 
intrinsic scatter matrices and the insensitivity of 
the fit parameter j3 in the vicinity of a c = 0. 

The blue (solid) curve is the result of restrict- 
ing the intrinsic scatter matrix to a range favored 
by the arguments in this section. Specifically we 
require the intrinsic scatter matrix color element 
to be within 3 standard deviations of the value 
determined by the fit to Equation (fT2|) . namely 
0.040 < <j c < 0.082. We also require <n < 1 
since the intrinsic scatter cannot exceed the mea- 
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sured width of the distribution, and £oc > since 
we expect a positive correlation when the color 
smearing includes a component that is uncorre- 
cted between the observed filter bands. While 
the exact values of these cuts are somewhat ar- 
bitrary, we have chosen a wide range to account 
for model uncertainties. Only 1650 of the 6380 
intrinsic scatter matrices satisfy these three crite- 
ria, but even with the restrictions a wide range 
of values are found. However, the median val- 
ues and the ranges that contains 68% of the data 
are a = 0.135±g;g?| and /3 = 3.19±oll- In addi- 
tion to the range of values that result from dif- 
ferent assumptions about the intrinsic scatter ma- 
trix, the statistical errors reported by SALT2mu 
are a a =0.011 and ap — 0.10. Adding these errors 
in quadrature, the final result is a = 0.1351°; ; i7 
and /3 = 3.261^24- F° r each entry in Figure [5] we 
compute the statistical probability that the entry 
is consistent with f3 > 4.1 using the computed sta- 
tistical error and assuming a normal distribution. 
The uncertainty in the intrinsic scatter matrix is 
computed by averaging the statistical probability 
for all the entries in Figure [5) From this calcu- 
lation, we conclude that the probability that our 
result is consistent with j3 > 4.1 is 2%. The prob- 
ability is higher than would be naively estimated 
from the value /3 = 3.26lo.24 because our calcula- 
tion includes the effect of the non-Gaussian tail of 
the distribution shown in Figure [SJd. 

6. Redshift Dependence of a and f3 in the 
SDSS data 

The parameters a and /3 are usually assumed 
to be constants, independent of redshift, but the 
possibility of a redshift dependence has been con- 
sidered previously by, for example, iKessler et al. 
(|2009bl ) . We might observe a redshift dependence 



of the fits for a and /3 for a number of reasons: 

• The properties of the SN, themselves, may 
evolve with redshift. 

• The model of Equation ((4]) may not be an 
exact model and the resulting fit parame- 
ters may be in tension between different SN 
populations. If the population is a function 
of redshift, the amount of tension may vary, 
resulting in different results for a and j3. 



• Our SN selection has biases, particularly at 
higher redshift where we tend to select the 
brighter SN. 

While we could directly fit for a(z) and fi(z), the 
likelihood that selection effects will play a signif- 
icant role suggests that a better approach might 
be to split the SDSS sample into 4 redshift bins. 
The results of the fits are shown in Figure [5] We 
have used the the same ensemble of intrinsic scat- 
ter matrices with u\ < 1, 0.040 < a c < 0.082, 
and £o c > that was used to construct Figured] 
The errors shown are the combined statistical and 
systematic error that arise from the uncertainty 
in the intrinsic scatter matrix. The results shown 
for a in Figure [6]^a) are consistent with the sam- 
ple average although there is an indication that 
the value of a may be rising as a function of red- 
shift. Figure [6^b) show a value of (5 that is 2.1cr 
lower in the lowest redshift bin, 2.2er higher in 
the range 0.22 < z < 0.32, and 3.7a lower in the 
highest redshift bin than the result for the com- 
plete sample. The accuracy in the determination 
of P as a function of redshift varies not only be- 
cause of the statistics but because highly reddened 
SN are too dim to meet our selection criteria at 
high redshift. The number of SN in each redshift 
bin and the number of highly reddened (c > 0.3) 
SN are: 36 and 4 (0.02 < z < 0.12), 149 and 2 
(0.12 < z < 0.22), 116 and (0.22 < z < 0.32), 
and 41 and (0.32 < z < 0.42), respectively. The 
reduced range of c at the higher redshifts reduces 
the accuracy of the fit and makes the fit consid- 
erably more sensitive to the amount of intrinsic 
scatter in the data. The highest-redshift sample 
consists of 41 SN with -0.319 < c < 0.126 that 
show very little correlation between the distance 
modulus and c, resulting in a low value of /3. The 
range of observed values of c is most likely primar- 
ily due to measurement error and intrinsic scatter 
resulting in the weak correlation between distance 
modulus and observed c. In principle, our formal- 
ism should recover the correct value of /3 when the 
correct measurement errors and intrinsic scatter- 
ing are included, but selection effects have signifi- 
cantly distorted the natural spread. In addition to 
the restricted color range, the absolute magnitudes 
of SN with 0.32 < z < 0.42 are 0.069 ±0.042 mag- 
nitudes brighter on average, assuming standard 
cosmological parameters (Q\ — 0.7 and w = — 1). 
As a consequence, our fit to the last redshift bin 
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should not be considered to be a good representa- 
tion of what would be obtained with an unbiased 
sample. A more robust result would be possible 
with an improved simulation of the selection ef- 
fects and augmenting the spectroscopically con- 
firme d sample with th ose identified photometri- 
cally (jSako et al.ll2011l ) to obtain a more complete 
sample. Selection effects should be much less im- 
portant for the midrange (0.22 < z < 0.32), where 
is j3 larger than the sample average by 2.2<r. While 
the lack of highly reddened SN (c > 0.3) degrades 
our ability to determine (3 in this redshift region, 
we can see if their absence biases f3 to higher val- 
ues by excluding SN with c > 0.3 in the low red- 
shift bins. The j3 values for the two lowest redshift 
regions change by less than 0.1 when the highly 
reddened SN are removed from the sample. 

7. Discussion 

The histograms in Figure[S]show the same trend 
that was observed in Table [TJ when the intrinsic 
scatter in color is large, the best fit value of (3 
will increase while a decreases slightly. The fact 
that our value for (3 is the highest value of those 
reported in $3] depends on our incorporation of a 
large intrinsic scatter in supernova color. We have 
not made any correction for the selection bias of 
our sample based on the results of the simulation, 
which is subject to the same selection criteria, and 
shows no significant bias in the determination of 
a and (3. 

The error on our result is dominated by an un- 
certainty in the intrinsic scatter matrix, but there 
are other potential problems that have not been 
taken into account. Certainly, our assumption 
that the intrinsic scatter is a constant, indepen- 
dent of any SN properties, is suspect given our 
lack of knowledge of the processes involved. In ad- 
dition, our model assumes that the data follow the 
relationship of Equation (j4]) with a fixed value of (3. 
However, we know that the v alue of Ry varies sig - 
nificantly in our own galaxy ( Schlaflv et alJuOlOl) , 
and we should expect to see at least some varia- 
tion in SN host galaxies. Variations in (3 produce 
effects similar to the intrinsic scatter in color, but 
the effect can be magnified in a magnitude lim- 
ited sample like SDSS where there are few highly 
reddened SN, which have higher weight in deter- 
mining (3. Variations in (3 could explain the small 



differences seen between the values of /3 in Figure 
[6] in the two lowest redshift bins. 

We have seen that neglecting intrinsic scatter in 
the S ALT2 c parameter results in a value of (3 that 
is biased low. The amount of bias depends on the 
amount of intrinsic scatter relative to the range in 
the variable c. Since the range in c decreases as 
the redshift increases, neglecting or underestimat- 
ing the amount of intrinsic scatter in c would be 
expected to result in a decreasing value of (3 as a 
function of redshift. The fits to the SDSS data 
as a function of redshift include intrinsic scatter 
in color so this effect should not be the explana- 
tion for the SDSS data — unless the scatter was 
underestimated. However, an overestimate of the 
intrinsic scatter in m x leads to similar biases and 
is more likely explanation for the lower value of f3 
that is observed in the highest redshift bin. 

A more problematic type of variation in (3 could 
arise from the color smearing. We have argued 
that color variations of SN must have some com- 
ponent of host galaxy extinction and also some 
intrinsic color variation, which could be due to 
variations in the explosion process or local envi- 
ronment or both. We have modeled the intrinsic 
color variation as color smearing with an effective 
value f3 — 0, but it is possible that each mechanism 
that contributes to the observed color requires a 
different correction to the distance modulus. Our 
model can only produce an effective value of f3 
that could be a function of color, redshift, host 
galaxy type or other parameters and could, there- 
fore, differ depending on the characteristics of the 
particular SN sample. 

Although we have not included correlations 
with galaxy types or spectral line features, the in- 
troduction of additional dependent variables into 
Equation (J4} is straightforward. In fact, having 
a more complete description of the supernova ex- 
plosion will not only provide a smaller dispersion 
for the numerator in Equation © , it will decrease 
the size of the intrinsic scatter that is required and 
decrease the importance of the uncertainty as to 
its form. 



Recently a number of papers (|Kim I 12011 ; 



Lago et al. 1 1201 It iMarch et al. 1 120111 ) have ex 



plored the technique of determining the intrin- 
sic scatter from the data. The approach of 



([March et al. 1 120111 ) is similar to ours in that it 



included our more general form for the intrinsic 
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error matrix but assumed that the off-diagonal 
terms involving Mb were zero. We have not con- 
sidered including the intrinsic scattering error as a 
fit parameter partly for simplicity and partly be- 
cause the underlying physics of the intrinsic scat- 
ter is still uncertain. However, the more rigorous 
statistical approach described by these authors is 
a possible direction for future work. 

We argue from the SDSS SN photometric data 
that there is likely some intrinsic scatter in SN 
color that is responsible for some of the scatter 
in the Hu bble diagram. A simi lar conclusion is 
reached by iFolev fc Kasenl (J2011I ) , based on much 
more detailed information involving differences in 
the line velocities measured from the SN spectra, 
which are modeled as arising from differences in 
viewing angle. While this type of detailed spec- 
tral information may not be available for all future 
SN surveys, an understanding of the underlying 
physical processes will be important for the most 
accurate treatment of the data. 

8. Summary 

We have described a new formalism to fit the 
parameters a and j3 that are used to determine 
the standard magnitudes of Type la supernovae. 
The determination of these parameters is accom- 
plished in a way that is independent of cosmology 
by introducing a few nuisance parameters to ac- 
commodate any reasonable cosmology. The for- 
malism also introduces a more general form for 
the intrinsic scatter which introduces a large un- 
certainty into the fit if the parameters describing 
the intrinsic scattering are not known accurately. 
We have shown by simulation that the mathemat- 
ical formalism is self-consistent, but have pointed 
out ways in which the model may fail to be a com- 
plete description of SN. 

We find that the SDSS data, when fit to 
the form of Equation (J4J are described by a = 

0.135±oloi? and t 



3.191q;24 by marginalizing 



over the uncertainty in the intrinsic scatter matrix. 
Our result relies on our conclusion that SNe are 
subject to a substantial color smearing in addition 
to reddening from host galaxy dust, as indicated 
by our fit to the color distribution where we deter- 
mined that <r c — 0.055 ± 0.007. The uncertainty 
in the parameters of the intrinsic scatter matrix 
results in a much larger error than is obtained if 



the parameters are assumed to be known. How- 
ever, even with the larger value of j3 and the larger 
error, we find that SDSS data differ at the 98% 
confidence level from (3 — 4.1, the value expected 
for extinction by dust in the Milky Way. 
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Fig. 4. — The distribution in SALT2 shape pa- 
rameter (xi) of 296 SDSS SN from our sample of 
343 SN that were also determined to be hosted 
by star-forming galaxies according to the analysis 
of ( Lampeitl et al.ll2010l ) is shown as a histogram 
and the dashed curve is Gaussian of width 1. The 
Gaussian is a good description of the core of the 
distribution but does not describe the tail towards 
low X\ . 
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Fig. 5.— The distribution of SALT2mu fit re- 
sults for a (a) and /3 (b) for a random distribu- 
tion of intrinsic scatter matrices that yield \ 2 P er 
degree of freedom equal to 1 is shown as the dot- 
ted red curve.. The solid blue curve shows the 
more restricted subset of error matrices that sat- 
isfy o-i < 1, 0.040 < er c < 0.082, and £ 0c > 0. The 
left axis shows the vertical scale for the red curve; 
the right axis is the scale for the blue curve. Of 
the 1650 entries in the restricted subset of intrinsic 
scatter matrices, there are 32 with a < 0, 15 with 
a > 0.4, 1 with j3 < and 15 with j3 > 5, which 
are not shown in the figures. The underflow and 
overflow entries result from cases where the en- 
tire distribution can be explained by the assumed 
intrinsic scatter of the data and the parameters 
a and j3 become indeterminate, resulting in large 
values and large estimated errors for the parame- 
ters. 
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Fig. 6.— The SALT2mu fit results for a (a) and /3 
(b) in 4 redshift bins using the same a distribution 
of intrinsic scatter matrices that yield chi-squared 
per degree of freedom equal to 1 and satisfy o\ < 
1, 0.040 < a c < 0.082, and £ 0c > 0. The results for 
the full sample are indicated by the dashed lines. 
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